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Evidence for pronounced quark loop effects in QCD 

Robert D. Mawhinney^* 

'"Department of Physics, Columbia University, New York, NY 10027, USA 

We have measured the hadron spectrum in lattice QCD, using staggered fermions, for (the quenched ap- 
proximation), 2 and 4 light degenerate dynamical quarks. In addition to earlier results involving extrapolations 
in valence quark masses for fixed dynamical mass, we also report results where we extrapolate in the dynamical 
mass for 4 flavors. We see a marked difference in the hadron spectrum for 2 and 4 flavors; the hadron spectrum 
is nearly parity doubled for 4 flavors, indicating smaller effects of chiral symmetry breaking. This pronounced 
effect in the hadron spectrum cannot be removed by a simple change in scale as the number of light quark flavors 
is changed. Further simulations at larger volume are needed to rule out finite volume effects. 



1. INTRODUCTION 

During the past two years, the lattice QCD 
group at Columbia has been exploring the effects 
of light (by today's standards) dynamical quarks 
on the hadron spectrum in zero temperature sim- 
ulations of QCD. Light quark effects are read- 
ily seen at finite temperature, where studies have 
shown the order of the QCD phase transition is 
sensitive to the number of light quarks, as ex- 
pected from theoretical arguments. Zero temper- 
ature simulations to date, including the ones re- 
ported here, have shown little difference between 
and 2 flavors of dynamical quarks. To see if this 
continues, we have done extensive simulations us- 
ing 4 flavors {Nf = 4). 

On our existing 256-node, 16 Gflop peak speed 
computer, we are essentially limited to lattices 
of size 16'^ X Nt for light quark mass simulations 
(ma 0.01, where m is the bare staggered quark 
mass). Using mp to set the scale gives us lattice 
spacings of ~ 2 GeV, for the lighter masses 
we simulate. (The exact value for the scale de- 
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pends on whether one uses nip extrapolated to 
zero quark mass or not.) This also makes our 
lattice about 1.5 fermi in the spatial directions. 
Quenched simulations show this introduces errors 
at the 5-10% level and similar systematic er- 
rors can be expected for the full QCD simula- 
tions, although complete studies await the arrival 
of more powerful computers. Since our volume 
is limited, we cannot extrapolate to the infinite 
volume limit and also not the continuum limit. 
Thus, our results are expected to be changed 
quantitatively as these limits are taken, but there 
is no compelling reason to expect a qualitative 
change. 

2. SIMULATIONS 

Table |l| lists the parameters for our simulations. 
It is important to note that the 2 flavor results 
employ an inexact algorithm Q , with errors pro- 
portional to (Ar)^. Since we have used the same 
step size, for a given dynamical mass, in both the 
2 and 4 flavor calculations and the 4 flavor cal- 
culations have a very high acceptance, the finite 
step size errors for the inexact algorithm are ap- 
parently very small. We have used a smaller con- 
jugate gradient stopping condition for the 4 flavor 
runs, since this enters into the accept/reject step, 
and we wanted to insure that this step is accurate. 

In these simulations we have measured hadron 
correlators for a variety of different source types 
and sizes with many different valence quark 
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Table 1 

Run parameters for the simulations whose results are reported here. The 2 flavor, ma = 0.015, 0.02 and 
0.025 results are from [3]. The longer thermalization time for the 0.025 run merely reflects measurements 
not being made from the beginning of the run. 





TOdynO = 0.01 


mdynO = 0.015 


mdynfl — 0.02 


TOdynO = 0.025 


A^^ = 2, /? = 5.7 


volume 


16^ X 40 


16^ X 32 


16^ X 32 


16^ X 32 


run length 


4870 


3010 


1425 


2830 


thermalization 


250 


250 


250 


530 


step size 


0.0078125 


0.0078125 


0.01 


0.01 


CG stopping condition 


1.01 X 10"^ 


1.13 X 10-=^ 


1.13 X 10-^ 


1.13 X lO-'^ 


iV/ = 4, /5 = 5.4 


volume 


16^ X 32 




16^ X 32 




run length 


4450 




2725 




thermalization 


250 




250 




step size 


0.0078125 




0.01 




CG stopping condition 


1.13 X 10-^ 




1.13 X 10~« 




acceptance 


0.95 




0.99 





masses (the masses which appear in the explicit 
Green functions for the hadron). This allows us 
to fit for excited states and check the stability 
of our effective mass plateaus [Q. We also have 
measured the chiral condensate for valence quark 
masses spanning 10 orders of magnitude which 
is useful as a probe of finite volume effects and 
chiral symmetry breaking. 

3. PARTIALLY QUENCHED RESULTS 

Results from some of these earlier simulations 
can be found in (5) and We studied the 
valence hadron spectrum for 0, 2 and 4 quark 
flavors with a single fixed dynamical quark mass 
mdynO = 0.01 for Nf ^2 and 4. For iV/ = 2 we 
used (3 = 5.7 and for Nf ^ A, [3 was 5.4. The 
parameters were chosen so that the rho mass, ex- 
trapolated to zero valence quark mass, was equal 
(we achieved equality at the 2% level) for 0, 2 
and 4 flavors. In choosing our parameters for the 
4 flavor simulation, we were helped by the earlier 
4 flavor work of . 

Figure |l] shows valence hadron masses for par- 
tially quenched 2 and 4 flavor simulations with 
rTidynfl = 0.01. The parity partners p, ai and iV, 



N' for 4 flavors show less splitting, particularly 
in the Wvaia — > limit, than for 2 flavors. 

The other major results of this work were (us- 
ing TUp to set the scale): 

1. TOjv and are very similar between and 
2 flavors. However, the mass of the parity 
partners, ttin'^ rriai , differs by 2-3a between 
these two cases, with ruN' and ruai closer 
to their parity partners for 2 flavors. 

2. For 4 flavors, mjsf/mp extrapolated to zero 
valence quark mass is 7% larger for 4 fla- 
vors than for 2 flavors (about 2a). The de- 
pendence of this quantitative difference on 
lattice volume and spacing deserves further 
study. 

3. Linear fits of to J to the valence quark mass 
do not go to zero as expected for a Gold- 
stone particle. This is consistent with the 
finite volume effects expected from a cutoff 
in the Dirac eigenvalue spectrum. In partic- 
ular, the mass of the Goldstone particle at 
zero quark masses is predicted to increase 
as the amount of chiral symmetry breaking 
decreases and this is seen in the data 101 . 
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Figure 1. Masses for the parity partners, p and 
oi, N and N' for a partially quenched calculation. 
Unfilled symbols are for 2 flavors, filled circles are 
for 4 flavors. All simulations had mdynO fixed at 
0.01. 



These results agree with the naive expecta- 
tion that adding more powers of the determinant 
(more quark fiavors) suppresses the small eigen- 
values of the Dirac operator, which decreases the 
effects of chiral symmetry breaking since the chi- 
ral condensate is given directly by the density of 
eigenvalues at the origin. The size of the effect, 
for the quark masses we are currently able to sim- 
ulate with, is certainly larger than anticipated, 
given the close agreement between the and 2 
flavor simulations. 



flavor simulations. 

Figure || shows the tt, p and nucleon masses 
for 2 and 4 flavors as a function of the dynamical 
quark mass. The unfilled symbols are the 2 flavor 
results; the filled symbols are for 4 flavors. The 
close agreement between the values for nip is due 
to our choice of parameters. The solid lines are 
fits to the 2 flavor results and show the extrap- 
olation to zero quark mass. The fits have large 
per degree of freedom (7, 3 and 5 for the tt, p 
and nucleon respectively). The reason for this is 
unknown (finite volume effects are a likely candi- 
date) and will hopefully be settled by future sim- 
ulations. The dashed lines extrapolate the two 
points we have for 4 flavors. They are consistent 
with the larger value for mj^jmp we found for 
4 flavors in the partially quenched hadron spec- 
trum. 




0.00 



0.01 



0.02 



4. FULL QCD RESULTS 

Having seen a marked decrease in the splitting 
between parity partners in the partially quenched 
hadron spectrum, we have been investigating the 
hadron spectrum in full QCD as a function oi Nf. 
To date, we have results for two dynamical mass 
values for 4 flavors and we can compare these to 
four dynamical mass values for 2 flavors. The 
three heavier masses for 2 flavors are results from 
m while the mdynfl = 0.01 is newer and includes 
3 times the data of the 0.01 point in [||. Table |^ 
gives some of the hadron masses for our 2 and 4 



Figure 2. The masses for the tt, p and N for 2 and 
4 flavors. The unflUed symbols are for 2 flavors 
and the flUed for 4. 



Figure || shows the p, oi, N, and N' masses 
for 2 and 4 flavors as a function of the dynamical 
quark mass. Again, unfllled symbols are for 2 fla- 
vors and flUed for 4. The 2 flavor fits have good 
values for the ai and N' . It is clear from the 
figure that the parity partners are not degenerate 
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Table 2 

Hadron masses for 2 and 4 flavor QCD. 





TOdynfl = 0.01 


mdynO — 0.015 


mdynO = 0.02 


mdynfl — 0.025 


iV/ = 2, /3 = 5.7 




0.249(2) 


0.293(2) 


0.349(2) 


0.388(1) 




0.435(5) 


0.455(8) 


0.501(7) 


0.551(4) 




0.564(13) 


0.594(22) 


0.638(18) 


0.755(32) 


iV 


0.663(8) 


0.685(10) 


0.781(10) 


0.839(6) 


iV' 


0.760(14) 


0.833(38) 


0.905(46) 


1.022(41) 


iV/ 4, /3 = 5.4 


TT 


0.292(5) 




0.357(3) 




P 


0.438(8) 




0.499(6) 




ai 


0.493(7) 




0.617(14) 




iV 


0.690(21) 




0.773(10) 




N' 


0.731(16) 




0.872(13) 





in the mdynO ^ limit. The dashed lines extrap- 
olate the 4 flavor results to zero quark mass and 
one sees the almost complete degeneracy in this 
limit. 




of whether these lattices are still in a conflning, 
chirally asymmetric phase. Figure ^ shows the 
staggered fermion chiral condensate, (xx), for 0, 
2 and 4 flavors as a function of the quark mass 
(the dynamical mass for 2 and 4 flavors). The 
solid lines are fits for and 2 flavors; the dashed 
line is an extrapolation for 4 flavors. The 4 fla- 
vor value at zero quark mass, 0.0023, is almost 4 
times smaller than the 2 flavor value, 0.00854(17). 
Using the Gell-Mann-Oakes-Renner formula for 
gives (using rrip at zero quark mass) 



U/mp{0) = 0.127 2 flavors 
U/mp{0) = 0.061 4 flavors 



(1) 
(2) 



The error on the 2 flavor result is 0.005, but this 
involves using fits of poor and hence the error 
is likely unreliable. However, the difference be- 
tween the 2 and 4 flavor results is large on the 
scale of the error. 



Figure 3. Masses for the parity partners, p and 
ai, N and N' . Unfilled symbols are for 2 flavors, 
flUed circles are for 4 flavors. 



This 



apparent degeneracy leads to the question 



5. CONCLUSIONS 

We have presented data showing a strong Nf 
dependence in the hadron spectrum for full QCD 
for fixed lattice volumes and spacings (in units of 
TOp). The results are consistent with the naive 
expectation that increasing the number of light 
quarks suppresses the part of the Dirac eigen- 
value spectrum responsible for chiral symmetry 
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Figure 4. The chiral condensate for 0, 2 and 4 fla- 
vors. For the quenched simulation, the horizontal 
axis is the valence quark mass. 



breaking. Clearly these simulations need to be 
extended to larger volumes and weaker couplings, 
which will require substantially more computer 
power. However, we are seeing a large efi^ect (al- 
most a factor of 2 change in /tt/ttIp) which is un- 
likely to be changed drastically by approaching 
the continuum limit. 

It is of some interest to ask if we have stum- 
bled into a new phase of QCD, possibly a lattice 
artifact. We have measured Wilson lines on our 
lattices, in the smaller 16 lattice spacing direc- 
tion, and find a small value consistent with a con- 
fined phase. We have no data to support a more 
exotic situation, confinement without chiral sym- 
metry breaking. A second question is the tem- 
perature of our lattices. For 4 flavors, estimating 
from existing data indicates the flnite tempera- 
ture phase transition should occur on lattices of 
size Nt = 10 - 12 for /3 = 5.4 §. Of course 
our lattices have A^^ = 16 and periodic boundary 
conditions in the spatial directions, but Ns = 16 
places this work in the confined phase. In addi- 
tion, the strong first order nature of the 4 flavor 
transition means that, even close to the critical 
coupling, the hybrid algorithm will be exploring 
the other phase extremely infrequently. 



It should be noted that extrapolating our two 
values for to^ for 4 flavors does not give rn^ — > 
as the dynamical quark mass goes to zero. We 
have seen this effect in the partially quenched sim- 
ulations, where it is consistent with finite volume 
effects. It deserves more study for the full simu- 
lations. In addition, it would be interesting to see 
if partially quenched Wilson spectrum measure- 
ments also showed the effects seen for staggered 
fermions. 

The author would like to thank Jim Sexton, 
Tony Kennedy and Robert Edwards for useful 
discussions and Catalin Malareanu for measuring 
the Wilson lines on these lattices. 
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